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Abstract 

We investigate an asteroseismic model of non-rotating paramagnetic neutron 
star with core-crust stratification of interior pervaded by homogeneous internal and 
dipolar external magnetic field, presuming that neutron degenerate Fermi-matter of 
the star core is in the state of Pauli's paramagnetic permanent magnetization caused 
by polarization of spin magnetic moments of neutrons along the axis of magnetic 
field of collapsed massive progenitor. The magnetic cohesion between metal-like 
crust and permanent-magnet-like core is considered as playing a main part in the 
dynamics of starquake. This process is thought of as impulsive release of magnetic 
field stresses, brought about by disruption of magnetic field lines on the core-crust 
interface resulting in the crust fracture, which is followed by tectonic shear displace- 
ments of crust against massive core. Focus is laid on the post-quake relaxation of 
the star by node-free torsional vibrations of highly conducting crustal solid-state 
plasma, composed of nuclei embedded in the degenerate Fermi-gas of relativistic 
electrons, about axis of magnetic field frozen in the immobile paramagnetic core. 
Two scenarios of these axisymmetric seismic vibrations are examined, in first of 
which these are considered as maintained by combined action of Lorentz magnetic 
and Hooke's elastic forces and in second one by solely Lorentz force. The damping 
is attributed to shear viscosity of crustal material. Based on the energy varia- 
tional method of magneto-solid-mechanical theory of elastic continuous medium, 
the spectral formulae for the frequency and lifetime of this toroidal mode are ob- 
tained and discussed in the context of theoretical treatment of recently discovered 
quasi-periodic oscillations of the X-ray outburst flux from SGR 1806-20 and SGR 
1900+14 as being produced by above seismic vibrations. 
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1 Introduction 



Studying seismic vibrations of neutron stars is among the most promising tools in search 
for dynamical laws governing continuum mechanics and macroscopic electrodynamics of 
extremely dense matter compressed by self-gravity to the nearly normal nuclear density 
at which microscopic composition of stellar material is dominated by neutron component. 
One of the most conspicuous features distinguishing this subclass of final stage (FS) stars 
of evolutionary track from their massive main-sequence (MS) progenitors is that they 
come into existence as sources of ultra powerful magnetic fields (Chanmugam 1992). It is 
hoped, therefore, that exploring vibrational behavior of spherical material masses capable 
of accommodating super strong magnetic fields we can gain some insight into the nature 
of the nuclear matter magnetism. 

In this connection it is appropriate to emphasize that contrary to liquid MS stars 
whose magnetic fields are generated in dynamo processes by flows of perfectly conducting 
flowing matter (energy supply of persistent convection in which comes from cental ther- 
monuclear reactive zones), in the FS stars, white dwarfs and pulsars, there are no nuclear 
energy sources to support persistent convection (e.g. Flowers & Ruderman 1977, Spruit 
2008). Such point of view invalidates applicability to FS stars of magneto-fluid- mechanical 
arguments and methods in the analysis of magnetic field effects such, for instance, as anti- 
dynamo theorems. Unlike the dynamo-generated fairly weak fields of liquid MS stars, the 
super strong and highly stable to spontaneous decay magnetic fields of isolated neutron 
stars (Bhattacharya & van den Heuvel 1991, van den Heuvel 1994, Bhattacharya 2002) 
are regarded as fossil (e.g. Ferrario and Wickramasinghe 2004). These are thought of as 
frozen in the immobile, non-convective, degenerate stellar matter possessing properties of 
elastically deformable solid. It is believed, therefore, that basic features of asteroseimol- 
ogy of strongly magnetized end-products of stellar evolution (white dwarf stars, neutron 
stars and still hypothetical quark stars) can properly be understood working from equa- 
tions of solid- mechanical theory of continuous media_j (McDermott et al 1988, Blaes et al 



x It may be appropriate to mention here investigations in nuclear physics on theoretical treatment 
of giant resonances in terms of vibrations of ultra fine piece of continuous nuclear matter in which it 
has been found that behavior of atomic nucleus in these fundamental modes of nuclear response bears 
strong resembles to spheroidal and torsional vibrations of elastically deformable solid sphere rather than 
drop of flowing liquid, (e.g. Bastrukov et al 2008b and references therein). One of the remarkable 
analytic inferences of these investigations is that the shear modulus of nuclear material (which is treated, 
therefore, as an elastic Fermi-solid) entering equations of nuclear solid-globe model, is found to be precisely 
coinciding with the pressure of degenerate Fermi-system of nucleons. This finding of nuclear physics 
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1989, Strohmayer 1991, Bastrukov 1996, Bastrukov et al 1999a, 1999b, Franco et al 2000, 
Bastrukov et al 2007 see, also, references therein). 

Many agree today that recent observations of quasi-periodic oscillations (QPOs) in the 
outburst X-ray flux of two soft gamma-ray repeaters, SGR 1806-20 and SGR 1900+14 
(Israel et al 2005, Watts & Strohmayer 2006) have offered a real opportunity of using 
the detected QPO frequencies for asteroseismological inferences about internal structure, 
mechanical and magnetic properties of material of these ultra magnetized neutron stars 
whose X-ray bursting activity exhibits many features in common with earth quakes (Chen 
et al 1996, Woods & Thompson 2006, Mereghetti 2008). The detected QPOs frequencies 
(measured in Hz), which are of particular interest for the subject of present work, are 



With all above in mind, in (Bastrukov et al 2007-2009) several models of post-quake vi- 
brational relaxation of a solid star have been investigated. Particular attention has been 
given to the regime of node-free or nodeless shear axisymmetric vibrations which remain 
less investigated in theoretical seismology as compared to the standing-wave regime of 
vibrations of both solid FS stars and such solid celestial objects as Earth-like planets (e.g. 
Lapwood & Usami 1981, Lay & Wallace 1995, Aki & Richards 2003). In particular, in 
(Bastrukov et al 2007b, 2008a), a case of the elastic-force-driven nodeless shear oscilla- 
tions, both torsional - t£ and spheroidal - S£, entrapped in the finite-depth crust AR of 
core-crust model of a solid star has been studied in some details with remarkable inference 
that dipole overtones of spheroidal and torsion nodeless elastic shear vibrations of crust 
against immobile core exhibit features generic to Goldstone soft modes. On the other 

suggests that neutron star can be thought of as a spherical mass of an elastic Fermi-solid of the normal 
nuclear density in which the energy of gravitational pull is counterbalanced by elastic energy stored in 
the incessant Fermi-motion of neutron quasi-particles in the potential of self-gravity (e.g. Bastrukov 
et al 2007a). Bearing this in mind and widely spread opinion that a neutron star can be regarded 
as an astrophysical counterpart of an atomic nucleus (e.g. Lipunov 1992, Bisnovatyi-Kogan 2002), it 
has been shown in (Bastrukov et al 1999a, 1999b) that equations of nuclear elastodynamics, originally 
introduced in the nuclear physics of giant resonances as fundamental equations of continuum mechanics 
of nuclear matter, can be successfully utilized for the study of neutron star pulsations. And in recent work 
(Bastrukov et al 2007a), it was shown that linearized equations of nuclear elastodynamics and canonical 
equations of linear theory of elasticity lead to identical results regarding frequencies of fundamental 
modes of node-free spheroidal and torsional vibrations of spherical mass of viscoelastic solid which are of 
particular interest for our present discussion. 



SGR 1900 + 14 : 28,54,84,155. 



SGR 1806 - 20 : 18, 26, 29, 92, 150, 625, 1840; 



(1) 
(2) 
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hand one can casts doubt on arguments of model presuming that detected QPOs have 
been brought about by vibrations restored by solely solid-mechanical force of Hooke's 
elastic stresses because such interpretation rests on the poorly justifiable assumption of 
the dynamically passive role of ultra strong magnetic field, as defining only direction of 
the axis about which the stellar matter undergo differentially rotational fluctuations, i.e., 
that magnetic field frozen in the star remains unaltered in this process. Bearing this 
in mind and assuming that the admixture of charged particles (basically electrons and 
protons) imparts to neutron-dominated Fermi-matter of these star mechanical properties 
of perfectly conducting material, it seems more plausible to expect that post-quake relax- 
ation should be strongly affected by Lorentz force of magnetic field stresses setting highly 
conducting material in rotational motion about axis of the dipole magnetic moment of 
the star. 

With this in mind, in (Bastrukov et al 2009a, 2009b), the above data have been 
analyzed in the model of a solid star undergoing global torsional nodeless vibrations 
about magnetic axis under the joint action of Lorentz force of magentic field stresses and 
Hooke's force of elastic shear stresses. And it was found that such a model provides 
fairly reasonable account of general trends in data on QPO frequencies for SGR 1900+14 
and data for SGR 1900+14 for frequencies from the range 30 < v < 200 Hz, but faces 
serious difficulties in interpreting the low-frequency vibrations with v — 18 and v = 26 
Hz in the SGR 1806-20 data. Also, the model of global torsional vibrations of solid 
star with froze-in homogeneous magnetic field leaves some uncertainties regarding the 
nature of vibrations with v = 625 and v = 1840 Hz. This last issue has been scrutinized 
in recent work (Bastrukov et al 2009b) from the standpoint of a solid star model with 
non-homogeneous poloidal magnetic field of well-known Ferraro's form and was found 
that these high-frequency QPOs can be properly interpreted as being produced by global 
Alfven node-free torsional vibrations of very high overtones. 

In the meantime, it has been argued some time ago in the studies of pulsar glitches 
as being produced by star quakes (Bastrukov et al 1996, 1997, 1999b) that post-glitch 
relaxation should most likely be dominated by Alfven oscillations of crustal electron- 
nuclear solid-state plasma (composed of immobile nuclei dispersed in degenerate Fermi- 
gas of relativistic electrons) about axis of magnetic field frozen in the immobile core 
(and, also, by less dense gaseous plasma expelled from the surface by starquake). From 
the view point of core-crust model of quaking neutron star (Franco et al 2000), there 
are two way of thinking why only peripheral region of the star can be set in Lorentz- 
force-driven Alfven vibrations. First belong to a case when strength of seismic tremor is 



4 



fairly week and perturbation sets in vibrations only finite-depth crustal region of the star, 
leaving massive core unaltered. The second way, when core material just incapable of 
sustaining Alfven vibrations. One of the main conditions for existing of such vibrations is 
the extremely large (effectively infinite) electrical conductivity of stellar substance (e.g., 
Chandrasekhar 1961, Mestel 1999), but this may not be the case for core material whose 
microscopic composition is dominated, according to canonical understanding of neutron 
star, by matter possessing properties of non-conducting degenerate Fermi-gas of non- 
relativistic neutrons. 

The physically meaningful framework for discussion of this second case provides the 
two-component, core-crust, model of quaking paramagnetic neutron star, pictured in 
Fig.l, that has been extensively discussed some time ago in works § (Bastrukov et al 
2001, 2002a, 2002b, 2003). In this model, the role of major source of super strong mag- 
netic field of the star is attributed to permanently magnetized massive core composed of 
degenerate neutron matter in the the state of Pauli's paramagnetic saturation in which 
magnetic moments of neutrons turn out aligned in the direction of frozen-in magnetic 
field that has been inherited by neutron star from massive progenitor and amplified in the 
implosive process of its gravitational collapse. Unlike the non-conducting core matter, the 
material of crust matter resembles metallic solid. This suggests that it is magnetic cohe- 
sion between metal-like crust and permanent-magnet-like core, mediated by lines of ultra 
strong magnetic field, should play main part in the the starquake dynamics. In this paper 
we adopt this latter point of view by considering starquake as impulsive release of mag- 
netic field energy by means of disruption of magnetic field lines on the core-crust interface 
and fracturing the crust by relived magnetic stresses. In some sense the magnetic field 
lines can be thought of as a super-hard piles endowing core-crust composition of neutron 
star with supplementary (to gravity forces) stiffness which is primarily responsible for its 

2 The main purpose of these works was to explore potential capabilities of equations of magneto-elastic 
dynamics (MED) suggested as fundamental equations of continuum mechanics and macroscopic electro- 
dynamics of non-conducting but highly magnetically polarized stellar matter that have been formulated 
in a manner of equations of magnetohydrodynamics underlying continuum-mechanical treatment of mo- 
tions of perfectly conducting liquid-state and solid-state plasmas pervaded by magnetic fields. It was 
found on the basis of MED equations that magnetically polarized paramagnetic neutron matter, with 
spin-magnetic moments of neutrons aligned along frozen-in magnetic field, can transmit perturbation by 
torsional magneto-elastic wave having for such a matter one and the same physical meaning as Alfven 
wave does for magneto-active liquid-state or solid-state plasmas. Also, it was shown that a neutron 
star whose material is in the state of paramagnetic permanent magnetization can sustain solely torsional 
oscillations driven by force of antisymmetric magneto-polarizational stresses. 
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PARAMAGNETIC NEUTRON STAR 




Figure 1: The internal constitution of two- component, core-crust, model of paramag- 
netic neutron star. The massive core is considered as a poorly conducting permanent 
magnet composed of degenerate Fermi-gas of non-relativistic neutrons in the state of 
Pauli's paramagnetic saturation caused by field-induced alignment of spin magnetic mo- 
ment of neutrons along the axis of uniform internal and dipolar external magnetic field 
frozen in the star on the stage of gravitational collapse of its MS progenitor. A highly con- 
ducting metal-like material of the neutron star crust, composed of nuclei embedded in the 
super dense degenerate Fermi-gas of relativistic electrons, is regarded as electron-nuclear 
solid-state magneto-active plasma capable of sustaining Alfven oscillations. 
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seismic stability to quake-induced tectonic shear displacements of crust against corqj. It 
is expected that quake-induced disruptions of magnetic field lines, along which the out- 
going radiation is developed, should be observed as glitches in the quiescent high-energy 
emission. These glitches are of substantially magnetic nature in the sense that they owe 
its origin to well-known effect of sudden, jump-like, increase in the magnetic field strength 
triggered by thermonuclear explosion. With this in mind we focus on reaction of the star 
on above quake-induced perturbation by node-free differentially rotational, torsional, vi- 
brations of crustal solid-state plasma about axis of magnetic field frozen in the immobile 
paramagnetic core. 

In section 2, a brief outline is given of the equations of solid-magnetics appropriate for 
the perfectly conducting viscoelastic continuous medium pervaded by a magnetic field. 
In section 3, the spectral formulae for the frequency and lifetime of this toroidal mode of 
non-radial nodeless vibrations driven by combined magnetic Lorentz and elastic Hooke's 
forces are obtained. The relevance of considered asteroseismic model to the detected 
QPOs in the X-ray flux during the flare of SGR 1806-20 and SGR 1900+14 is assessed in 
section 4. The newly obtained results are briefly summarized in section 5. 



3 The internal constitution of paramagnetic neutron star model under consideration does not exclude 
the existence of a fairly thin inner crust with average density several times less than the normal nuclear 
density (neutron drip-line) at which micro-composition of matter is dominated, as we know from nuclear 
physics, by paired neutrons. It must be stated clear, however, that such quasi-boson matter composed of 
paired neutrons has nothing to do with fairly dilute and flowing quasi-boson matter like supcrfluid Fermi- 
liquid composed of 3 He atoms. According to current theory, the emergence of superfluidity in flowing 
boson material substances belongs to phase transition of macroscopic fraction of bosons to the state of 
Bose- Einstein condensation (BEC). One of the most striking features of such phase transition is that at 
the temperature of condensation the internal pressure of BEC state vanishes. The fact that Bose-matter 
composed of paired neutrons in BEC state is unable to withstand the pressure of self-gravity (contrary to 
the degenerate neutron Fermi-matter whose zero-temperature pressure of degeneracy counterbalancing 
the self-gravity pressure is central to the very notion of neutron star) suggests that such a matter, if exists, 
can only insufficiently contribute to the total mass budget of neutron star, otherwise the very notion of 
this compact object is loosing its original meaning. In the starquake of paramagnetic neutron star model 
under consideration, this inner crust material is thought of as operating like a lubricant facilitating non- 
compressional differentially rotational shear displacements of crust relative to much denser matter of 
massive core. 
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2 Governing equations of solid-magnetics 



We study an asteroseismic model of non-rotating neutron star responding to quake- 
induced perturbations by shear (substantially non-compressional, bp = —pV k u k = 0) 
differentially rotational fluctuations of material displacements Ui, locked in the finite- 
depth crust of the density p. And it is believed that dynamics of this vibrational process 
can be properly modeled by equations of classical solid mechanics 

pUi = VfcTjfc + VfcCTifc + Vfc7T ifc , VfcWfe = (3) 

presuming that Hooke's elastic stresses a ik and Newton's viscous stresses ir ik are described 
by linear constitutive equations of the form 

a ik = 2pu ik , u ik = ^[ViU k + V k Ui\, (4) 

n ik = 2r]u ik , u ik = -[VjUfc + V k Ui] (5) 

where p stands for the shear modulus, 77 for shear viscosity and u ik is the tensor of shear 
strains or deformations. The central to our further discussion is the tensor of fluctuating 
magnetic field stresses 

r tk = ]- [Bi5B k + B k SBi - BjSBjS^} , (6) 

47T 

5Bi = V fc [ui B k - u k Bi] (7) 

As in our previous work (Bastrukov et al 2009a), we consider model with homogeneous 
internal magnetic field whose components in spherical polar coordinates read 

B r = B cos 9, B e = -B sin 9, B^ = (8) 

and external dipolar magnetic field is described by B = V x A, where A = [0, 0, = 
m s /r 2 ] is the vector potential with the standard parametrization of the dipole magnetic 
moment m s = (1/2)1? i? 3 of star of radius R and by B is understood the magnetic field 
intensity at its magnetic poles, B = B p . 

2.1 The energy method 

This method of computing frequency of shear vibrations rests on the equation of energy 
balance 



d f pit 
dt 



/pil 2 f 
~Y dV = ~ J ^ Tik + aik + ^ ^ dV ^ 



which is obtained by scalar multiplication of equation of magneto-solid-mechanics ([3]) with 
Ui and integration over the volume of seismogenic layer. From the technical point of view, 
the shear character of material distortions brought about by forces under consideration 
owes its origin to the symmetric form of stress-tensors in terms of which these forces are 
expressed. The key idea of this method consists in using of the following separable form 
of material displacements 

Ui(r,t) = fli(r)a(t) (10) 

where aj(r) is the time- independent solenoidal field and amplitude a(t) carries information 
about temporal evolution of fluctuations. Thanks to this latter substitution, the all tensors 
of fluctuating stresses and strains take similar separable form 

T ik (r,t) = [fifc(r) - ^rjj(r)8ik]a(t), (11) 

f*(r) = -jj-[Bi(r) 6 fc (r) + B k {r), 6;(r)], (12) 

bi(r) = V k [a<(r) B k (r) - a k (r) B^r)], (13) 

cr ifc (r) = 2fia ik (r)a(t), TT ik (r) = 2rj a ik (r)a(t), (14) 

a»fc(r) = -[V»ajfc(r) + V^a^r)]. (15) 

On inserting (fTU|) - (|T4"]) in the integral equation of energy balance (jUj) we arrive at equation 
for a(t) having the well-familiar form 

d£ „ n Ma 2 Ka 2 ^ Va 2 

- = -2T, s = — + — r=— (16) 

Ma + Va + JCa = 0, (17) 
a(t) = a exp(— tjr) cos(Ot), (18) 
fi 2 = W 2 



2 K 2M 
w = A?' T " • (19) 



where the inertia Ai, viscous friction T> and stiffness K, of oscillator are given by 

M = J P (T)ai(T)ai(T)dV, /C = /C e + /C m , (20) 

K e = 2 J fi(r) a ik {r) a ik {r) dV, (21) 

lC m = J f ik (r) a ik (r) dV, (22) 

V = 2 |r/(r)a iA; (r)a 4fc (r)dV. (23) 

The physical significance and practical usefulness of the outlined energy method is that it 
can be efficiently utilized in the study of non-radial seismic vibrations of more wide class 
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of solid degenerate stars like white dwarfs stars (Lou 1995, Molodtsova et al 2010) and 
ultra dense quark-matter stargj whose material, as is expected (Xu 2009 and references 
therein), can too be thought of as an elastic solid extremely robust to compressional 
distortions. In (Bisnovatyi-Kogan 2008, Heyvaerts et al 2009) it is proposed that ultra 
powerful luminosity of currently monitored 7-ray sources could be related to torsional 
vibrations of neutron and quark stars. 



2.2 Material displacements in the crust undergoing node-free 
torsional vibrations 

In this subsection we pause on detailed derivation of material displacements in the crust 
undergoing non-radial torsional oscillations in the node-free regime. In this mode, the 
velocity of fluctuating material flow (the rate of material displacements) is described by 
general formula of rotational motions 

*v(r,t) = u(r,t) = [n(r,t)xr], (24) 
n(r, t) = [V x <fv(r, t)] = [V x u(r, £)] = $(r, t) (25) 

However, unlike a case of rigid-body rotation, in which the angular velocity is a constant 
vector, in a solid mass undergoing axisymmetric differentially rotational vibrations the 
angular velocity fi(r, t) is the vector-function of position (Bastrukov et al 1999a, 2007a) 
which can be represented as 

n(r, t) = $(r, t) = 0(r) a(t), <f>(r) = [V x a(r)] (26) 

In the regime of node-free vibrations in question, a(r) is described by the divergence- 
free odd-parity, axial, toroidal field which is one of two harmonic solenoidal fields of 
fundamental basis (Chandrasekhar 1961) obeying the vector Laplace equation V 2 a = 0. 
This field can be expressed in terms of general solution of the scalar Laplace equation as 
follows 

a(r) = a t (r) = V x [r X (r)] = [V*(r) x r] (27) 

V 2 x(r) = 0, (28) 

x(r) = [A e r e + B e r- £ - 1 ]P e (C), ( = cosfl (29) 

4 The average density of quark star matter is expected to be comparable with the average density of 
nucleon p « 4.3 x 10 17 g cm -3 , that is, three order of magnitude greater than average density of normal 
nuclear matter p sw 2.8 x 10 14 g cm~ 3 . 
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where Pe(() is the Legendre polynomial of multipole degree i. It follows that the angular 
field 4>(r) is the poloidal vector field 



0(r) = [V x a t (r)] = V x V x [rx(r)] (30) 

= V^Z + ^r-^P^C), (31) 

Ai = A t (£ + l), B e = B e l (32) 

And this field is irrotational, V x 0(r) = 0. 



Figure 2: Material displacements of crustal matter about the dipole magnetic moment 
axis of paramagnetic neutron star undergoing nodeless differentially rotational, torsional, 
vibrations in quadrupole and octupole overtones. 



As was stated, we study a model of differentially rotational vibrations of peripheral 
finite-depth crust against immobile core. In this case, the arbitrary constants At and Be 
can be uniquely eliminated from two boundary conditions: (i) on the core-crust interface 
r = R c 



and (ii) on the star surface r = R 



U<f>\ r =R = [<& X Ry r=i? , 

* = a(t)V fl P/(C), K = e r R 



where 



d 1 d 



(33) 



(34) 
(35) 



(36) 



09 9 sm6d(j)_ 

The no-slip condition on the core-crust interface, r = R c , reflects the fact that the am- 
plitude of differentially rotational oscillations gradually decreases down to the star center 
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and turns into zero on the core. The boundary condition on the star surface, r = R, 
is dictated by symmetry of the general toroidal solution of the vector Laplace equation 
which then is tested to reproduce the moment of inertia of a rigidly rotating solid star 
(Bastrukov et al 2008a). The above boundary conditions lead to the coupled algebraic 
equations 

A t R[7 x + B £ R^- 2 = 0, A t R l + BiR-t- 1 = R (37) 

whose solutions are 

At = Aft, B t = -NtRT\ ^ = R2l+l _ Rf+1 - (38) 

In spherical polar coordinates, the nodeless toroidal field has only one non-zero azimuthal 
component 



a r = 0, ag = 0, (39) 
J , *i 1 (1 _ c2) i/ 2 ^0 (4Q) 



Air 1 



The snapshot of material node-free displacements in the crust undergoing torsional oscilla- 
tions against immobile core of paramagnetic neutron star under consideration is pictured 
in Fig.2 for quadrupole, £ = 2, and octupole £ = 3, overtones of this axial mode. 

3 Spectral formulae for the frequency and lifetime of 
seismic vibrations 

The computation of integrals defining mass parameter Ai, parameter of vibrational rigid- 
ity /C and viscous friction T> [which has been presented in some details elsewhere (Bas- 
trukov et all 1999, 2007a, 2009)] are quite lengthy but straightforward and, therefore, are 
not presented here. The mass parameter can be conveniently represented in the form 

M - i ^ R5 wfm^) m{e) ' (41) 

m(£) = (1 - A 2m )~ 2 (42) 



x 



l-(2£ + 3)\ 2e+1 + 



(2£ + l) 2 2e+3 1i + 3 > 2 (2M-i) 

+ ^7^i ~27^I A 

\ = ^ = l-h, h=^, (43) 
AR = R-R C , 0<A<1. (44) 
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The A-terms in the above and foregoing equations emerge as result of integration along 
the radial coordinate from radius of the core-crust interface r = R c to the star radius, 
r = R. The integral coefficient of viscous friction is given by 

V = 4n V R 3i -^^d(£) } (45) 

d(£) = (1 - A 2£+1 )- 



(* + 2) 
(£-1) 



A 



2£+l 



(46) 




Figure 3: The fractional lifetime of torsion nodeless oscillations of the neutron star crust 
damped by force of viscous shear stresses as a function of multipole degree £ computed 
at indicated values of the fractional depth h of peripheral seismogenic layer. 

For the lifetime we obtain 

/ \ 2r„ m{£) R 2 7] 

tM= (2£ + 3)(£-l)W)> T » = -> V= ~p (47) 

In Fig.3, the fractional lifetime is plotted as a function of multipole degree £ with indicated 
values of fractional depths h = AR/R. It shows that the higher £ the shorter lifetime. 
It is easy to see that in the limit, A = (R c /R) — > 0, we regain the spectral formula for 
lifetime of global torsional nodeless vibrations of solid star (Bastrukov et al 2003, 2007a) 

2tV R 2 _ fj 

TM= (2£ + 3)(£-lV ^ = ^' V = -p (48) 

in which by f v is understood, in this latter case, the average kinematic viscosity of the star 
matter as a whole; the extensive discussion of this transport coefficient can be found in 
(Shternin 2008). For the node- free torsional oscillations of solid star, the last equation has 
one and the same physical significance as the well-known Lamb formula does for the time 
of viscous damping of spheroidal node-free vibrations which in the context of neutron 
star pulsations has been extensively discussed some time ago by Cutler and Lindblom 
(1987). Regarding the problem under consideration we cannot see, however, how the 
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obtained formulae can be applied to observational data on QPOs in SGRs. Nonetheless, 
their practical usefulness is that they can be utilized in the study of a more wide class of 
celestial objects such as Earth-like planets (Lay & Wallace 1995, Aki & Richards 2003) 
and white dwarf stars. 

From above it is clear that the integral coefficient of elastic rigidity /C e of torsional 
vibrations has analytic form similar to that for coefficient of viscous friction £>, namely 

0(02 ~ 1) 



2£+l 



k e (£), 



1-A 



2£+l\-l 



+ 2) , 2 e+i 



-A 



(49) 
(50) 



(£-iy 

The frequency as a function of multipole degree £ of node-free elastic vibrations in question 
v e (£) (measured in Hz and related to angular frequency as uj e (£) = 2ixv e (£) = K e j M) is 
given by 



u 2 (£) = u 2 [(2£ + 3)(£~l)] 1 / 2 



k e (£) 



2-RV e = — , 

R' 



h, h 



AR 
~R~ 



(51) 
(52) 



It is easy to see that in the limit A — > 0, we regain spectral formula for the frequency 
of global torsional oscillations, having the form of equation (IBTj) with k e (£) = m(£) = 1. 
Understandably that in this latter case all material characteristics belong to the star as 
a whole. 

The magneto-mechanical stiffness of Alfven vibrations K, m can conveniently be written 

as 

£{£ 2 -1){£+1) 
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This leads to the following two-parametric spectral formula 
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(55) 

v^p (56) 

In the forward asteroseismic analysis of QPO data relying with on this latter spectral 
formula, the Alfven frequency, z/4 and the fractional depth of seismogenic zone, h, are 
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regarded as free parameters which are adjusted so as to reproduce general trends in the 
observed QPOs frequencies. 

Alfven toroidal mode 
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Figure 4: Fractional frequency and period of nodeless torsional magneto-solid-mechanical 
oscillations, toroidal Alfven mode - o<4> entrapped in the neutron star crust as functions 
of multipole degree I computed at indicated values of the fractional depth h of peripheral 
seismogenic layer. The value h = 1 corresponds to global torsional oscillations excited in 
the entire volume of the star. Here va = uja/2tt, where uja = va/R with va = B/y/Anp 
being the velocity of Alfven wave in crustal matter of density p and Pa = 2ti/ua- 

In Fig.4, the fractional frequencies and periods of this toroidal Alfven mode as func- 
tions of multipole degree £ are plotted with indicated values of fractional depth of the seis- 
mogenic layer h. Remarkably, the lowest overtone of global oscillations is of quadrupole 
degree, I — 2, whereas for vibrations locked in the crust, the lowest overtone is of dipole 
degree, i = 1, as is clearly seen in Fig.5. This suggests that dipole vibration can be 
thought of as Goldstone's soft mode whose most conspicuous property is that the mode 
disappears (the frequency tends to zero) when key parameter regulating the depth of 
seismogenic zone A — > 0. 
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Figure 5: Fractional frequency of nodeless torsional Alfven oscillations of indicated over- 
tones £ as a function of the fractional depth h of peripheral seismogenic layer. The 
vanishing of dipole overtone in the limit of h — > 1, the case when entire mass of neutron 
star sets in torsional oscillations, suggests that dipole vibration possesses property typical 
to Goldstone's soft modes . 

4 Application to QPOs in the outburst X-ray flux 
from SGR 1806-20 and SGR 1900+14 

The basic physics underlying current understanding of interconnection between quasi- 
periodic oscillations of detected electromagnetic flux and vibrations of neutron star has 
been recognized long ago (e.g. Boriakoff 1976, van Horn 1980, Blaes et al 1988). Owing 
to the effect of strong flow-field coupling, which is central to the propagation of Alfven 
waves, the quake induced perturbation excites coupled vibrations of perfectly conducting 
solid-state plasma of the crust (as well as gaseous plasma of magnetar corona expelled 
from the surface by outburst) and frozen-in lines of magnetic field. Outside the star the 
vibrations of magnetic field lines are coupled with oscillations of gas-dust plasma expelled 
from the star surface by quake. And it is these fluctuations of outer lines of magnetic field, 
operating like transmitters of beams of charged particles producing coherent (curvature 
and/or synchrotron) high-energy radiation, are detected as QPOs of light curves of the 
SGRs giant flares. 

In applying the obtained spectral formulae to the frequencies of detected QPOs we 
examine two scenarios, namely, when quake-induced torsional vibrations are restored by 
joint action of Lorentz magnetic and Hooke's elastic forces and when oscillations are of 
pure Alfven's nature, that is, produced by torsional seismic vibrations of crust against 
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core under the action of solely one Lorentz force of magnetic field stresses. 

4.1 Crust vibrations driven by combined action of Loternz mag- 
netic and Hooke's elastic forces 

In this case, the asteroseismic analysis of detected QPOs rests on the three-parametric 
spectral formula 

v\t)\v A , v e , h] = u 2 m (£)[u A , h] + vl{t)\v e , h] (57) 

The suggested theoretical £-pole specification of the detected frequencies is presented in 
Fig. 6 for SGR 1900+14 and in Fig. 7 and Fig. 8, exhibiting remarkable correlation 
between depth of seismogenic zone the fundamental frequencies of coherent magnetic and 
elastic oscillations - the larger AR, the higher basic frequencies of Alfvenic v& and elastic 
v e vibrations. It is seen from computations for SGR 1806-20, that reasonable fit of data 
can be attained with h = 0.2 (for the star model with radius 20 km, AR = 2 km) and 
with h = 0.4 (AR = 5 km). 
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Figure 6: Theoretical fit of the QPOs frequency in the X-ray flux from SGR 1900+14 
on the basis of three-parametric theoretical spectrum of frequency of torsional seismic 
vibrations in the crustal region of indicated fractional depth. 

It is worth emphasizing that at above values of h, the obtained here tree-parametric 
spectral formula much better match the data as compared to that for global, in the entire 
volume, vibrations studied in (Bastrukov et al 2009a). On this ground we conclude, if the 
detected QPOs are produced by seismic vibrations of peripheral region of the star under 
coherent action of Lorentz and Hooke's forces, then the depth of seismogenic layer AR 
should be quite large, somewhere in the range 0.2R < AR < OAR. 




5 10 15 20 25 30 35 
I 



17 



200 




100 150 200 250 300 350 400 
I 



Figure 7: The same as Fig.6, but for SGRs 1806-20 with h=0.2. 

4.2 Seismic vibrations of peripheral layer under the action of 
Loternz magnetic force 

In this case £-pole specification of detected QPOs frequencies rests on the two-parametric 
spectral formula for frequency of purely Alfven toroidal mode 

v 2 (t) = v 2 j£[v A ,h\). (58) 

The results presented in Fig.9 and Fig. 10 show that at indicated input parameters, 
i.e., the Alfven frequency Va and the fractional depth of seismogenic layer h = AR/R, the 
model too adequately reproduces general trends in the data with fairly reasonable £-pole 
specification of overtones pointed out by integer numbers along x-axis. It is seen that the 
low-frequency QPOs in data for SGR 1806-20, are interpreted as dipole and quadrupole 
overtones: ^(o^i) = 18 and £(00.2) = 26 Hz. And the high-frequency kilohertz vibra- 
tions with 627 Hz and 1870 Hz are unambiguously specified as high-multipole overtones: 
v(oa\ =i2 ) = 627 Hz and ^(0^=122) = 1870 Hz. 

However, in this latter scenario of Lorentz-force-dominated vibrations the best fit of 
data is attained at fairly large value of fractional depth, h = 0.6, which is much larger than 
the expected depth of the crust. In our opinion, this result may be regarded as indication 
to that the detected QPOs are formed by coherent vibrations of crustal solid-state plasma 
and plasma of magnetar corona. 
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Figure 8: The same as Fig.6, but for SGRs 1806-20 with h=0.4. 

5 Concluding remarks 

The neutron stars, both pulsars and magnetars, are seismically active compact objects 
which are formed, as is widely believed, in the magnetic-flux-conserving supernova collapse 
of the cores of massive stars coming into existence as sources of super strong magnetic 
field with a fairly complex internal constitution. The internal micro-composition of these 
compact objects has been and still remains among the most active investigations on 
search for adequate equations of state of neutron star matter (e.g. Weber 1999, Lattimer 

6 Prakash 2007, Sedrakian 2007, Haensel, Potekhin & Yakovlev 2007). 

Together with this development, the past three decades have seen increasing recogni- 
tion that important insight into interior of neutron stars and properties of super dense and 
ultra magnetized matter provides asteroseimology which seeks to explain fast variability 
of their luminosity as being produced by seismic vibrations. Adhering to this attitude, we 
have investigated seismic response of core-crust model of quaking paramagnetic neutron 
star, pictured in Fig.l, by node- free torsional crust-against-core vibrations. In this model 
the source of ultra powerful and highly stable to spontaneous decay magnetic field is asso- 
ciated with massive core, which is thought of as an ultra powerful spherical magnet. The 
physical nature of seismic stability of core-crust coupling is attributed to magnetic cohe- 
sion between core and metallic like crusto The fact that impulsive release of magnetic 



5 At this point it may be appropriate to note that considered model permits co-existence of two sources 
and two compound parts of neutron star magnetism (current-carrying flows in the crust and field-induced 
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Figure 9: Theoretical description (lines) of detected QPO frequencies (symbols) in the 
X-ray flux during the flare of SGRs 1806-20 and 1900+14 as overtones of pure Alfven 
torsional nodeless oscillations of crustal magneto-active plasma under the action of solely 
Lorentz restoring force. 

5000 
1000 

N 

X 

> 100 



"0 50 100 150 

I 

Figure 10: Same as Fig.9 but for SGRs 1806-20. 

energy in the starquake, associated with disruption of magnetic field lines on the core- 
crust interface followed by fracturing the crust by relived magnetic stresses, suggests that 
model under consideration is in line with current understanding of X-ray bursting activity 
of soft gamma-ray repeaters. With this in mind, the obtained spectral formulae have been 
applied to asteroseismic analysis of data on quasi-periodic oscillations of the outburst X- 

spin paramagnetism of core) and consistent with the idea about strengthening of magnetic field intensity 
in young neutron stars, like magnetars, due to thermoelectric effect in the crust (Urpin & Yakovlev 1980, 
Blandford et al 1983). And also with the idea that damping of crustal currents, thought of as second 
source of magnetic field of pulsars and magnetars, serves as one of the main factors responsible for the 
decay of magnetic field of neutron stars (e.g. Goldreich & Reisenegger 1992, Pons & Geppert 2008). 
It is expected that magnetic field decay should manifest itself in the lengthening of periods of Alfven 
oscillations. 
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ray emission from SGR 1900+14 and SGR 1806-20. In so doing we have investigated 
two scenarios of post-quake relaxation of paramagnetic neutron star by torsional vibra- 
tions of crust against immobile massive core distinguishing restoring forces of vibrations. 
In first case, the analysis of data on QPOs frequencies has been based on assumption 
that detected QPOs owe their existence to torsional nodeless vibrations restored by joint 
action of Lorentz magnetic and Hooke's elastic forces and we have found that obtained 
three-parametric spectral formula provides much better fit of data than two-parametric 
frequency spectrum of the early considered model of global magneto-elastic vibrations 
(Bastrukov et al 2009a). In second case of vibrations restored by solely Lorentz force of 
magnetic field stresses we found that obtained two-parametric frequency spectrum can too 
be fairly reasonably reconciled with detected QPOs frequencies. With all above in mind 
we conclude that main part in quake-induced torsional vibrations of crustal solid-state 
plasma about axis of magnetic field frozen in the immobile core plays Lorentz restoring 
force of magnetic field stresses and that observed frequencies of QPOs can be interpreted 
as an indirect evidence for permanent magnetization of super dense matter in the deep 
interior of magnetars. 

The authors are grateful to Dima Podgainy (JINR, Dubna) for helpful assistance, 
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